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Abstract 

We show that quasi-long-range (QLR) pair correlations that decay asymptotically with scaling 
r —(d+i) j n ^-dimensional Euclidean space M. d , trademarks of certain quantum systems and cosmo- 
logical structures, are a universal signature of maximally random jammed (MRJ) hard-particle 
packings. We introduce a novel hyperuniformity descriptor in MRJ packings by studying local- 
volume- fraction fluctuations and show that infinite- wavelength fluctuations vanish even for packings 
with size- and shape-distributions. Special void statistics induce hyperuniformity and QLR pair 
correlations. 



Quasi-long-range (QLR) pair correlations with asymptotic scaling r 



-(<*+!) 



are unique 



trademarks of noninteracting spin-polarized fermionic ground states ground-state 
liquid helium {3], and the Harrison-Zeldovich spectrum of the early Universe {4J. Such cor- 
relations also arise in maximally random jammed (MRJ) packings of monodisperse spheres 



j^j, which are prototypical glasses [6| possessing maximal disorder among all jammed pack- 
ings with diverging elastic moduli [7H9(. These packings were once thought to describe the 
structure of liquids [10j, which typically possess pair correlations decaying exponentially 
fast . It is thus unclear why special QLR correlations occur in MRJ packings and 

whether this behavior extends to packings of particles with size- and shape-distributions, 



which possess nontrivial material properties 14, Il5j. 



Noting that MRJ packings are structurally rigid with a well-defined contact network, 
Torquato and Stillinger conjectured [?| that all strictly jammed, i.e., globally incompressible 
and nonshearable, saturated packings of monodisperse spheres in d-dimensional real space 
M. d are hyperuniform, i.e., infinite- wavelength local- number-density fluctuations vanish 
a proposition for which no counterexample has been found to date. Although this conjec- 
ture made no mention of QLR pair correlations, further work on 3D MRJ monodisperse 
sphere packings js| established that the structure factor S(k) approaches zero linearly as 
the wavenumber k — > 0, inducing a QLR power-law tail r -4 in the pair correlation function 
92{ r )i which is proportional to the probability density of finding a particle within a radial 
distance r from a reference particle. However, neither study provided a quantitative ex- 
planation for the presence of QLR correlations in these systems, and the analysis does not 
blindly extend to MRJ packings of nonspherical and/or polydisperse particles in which the 
shape- and size-distributions of the particles are crucial. Also, it is not obvious why MRJ 
packings are hyperuniform since they lack the crystalline long-range order and Bragg peaks 
of periodic structures 16[ and yet exhibit QLR pair correlations. These unusual character- 
istics have important implications for understanding the structural properties of MRJ and 
glassy material systems 17|, LL8| . We note that hyperuniform disordered point patterns have 
been used to create new materials with unusual optical properties [19[ and large, complete 



photonic band gaps 



20|. 



In this Letter, we characterize universal properties of polydisperse and/or nonspherical 
MRJ hard-particle packings. Specifically, we provide strong support for the claim that all 
saturated, strictly jammed particle packings, which may be polydisperse and/or nonspher- 



ical, possess vanishing infinite-wavelength local-volume-fraction fluctuations, a recently 



introduced extended notion of hyperuniformity [11] . We thus provide the appropriate novel 
framework and descriptor to study universal QLR pair correlations in general MRJ pack- 
ings. By examining binary circular disk, ellipse, and superdisk packings in IR 2 , we show that 
infinite-wavelength density fluctuations associated with the particle centers do not vanish 
due to local microstructural inhomogeneities. However, infinite-wavelength local-volume- 
fraction fluctuations do vanish, so that the appropriate structural descriptor of MRJ pack- 
ings is the two-point probability function S2 (r), defined as the probability of finding both of 
two arbitrary points either within particles or in the exterior void space. Importantly, our 
methodology contains all previously-reported results for monodisperse MRJ sphere packings 
5|, |7J] as a special case, meaning that we provide a completely general means of understand- 
ing QLR correlations in jammed hard-particle packings. We show that the rigidity of the 
contact network places strong constraints on the size- and spatial-distributions of the local 
voids, leading to hyperuniform QLR pair correlations that compete with the maximal ran- 
domness of the packings. This observation supports the important notion that the topology 
and geometry of the void space 



2l| are more basic than the particle space [ijj]. 
It is sufficient for our purposes to consider 2D packings of binary circular disks, ellipses 
and superdisks, but our arguments concerning the void space are general enough to incorpo- 



rate higher-dimensional packings. We use the Donev-Torquato-Stillinger algorithm 22j to 
generate jammed packings of hard particles in IR 2 . Particles of two different sizes undergo 
event-driven molecular dynamics with periodic boundary conditions while simultaneously 
growing at a specified rate and fixed size ratio, resulting in a final configuration that is 
essentially strictly jammed with a packing fraction <fi « 0.8475. The configurations that 
we study contain a small concentration (~ 2.5%) of "rattler" particles, which are free to 
move within some small cage. Saturation of the packing places an upper bound on the cage 



size, and the rattler locations are roughly uniformly distributed in the simulation 
moving these rattlers breaks hyperuniformity by lifting the saturation constraint 



dox. Re- 



1 



these 



particles are thus kept in our configurations, which, by numerical protocol design, are close 
approximants to the MRJ state 5j. We have chosen particle concentrations 7 sm aii = 0.75 
and 7iarge = 0.25 with size ratio (3 = 1.4 and have found that varying these parameters does 
not affect hyperuniformity. 

Hyperuniformity is determined by pair separation distances within either a point pattern 



or a heterogeneous medium. Analysis of point patterns thus requires the pair correlation 
function g 2 or the associated structure factor S(k) = 1 +p${g2(r) — l}(k) with $ the Fourier 



transform. Using these 
variance o~ 2 N {R) in M 2 as 



unctions, one can write the asymptotic behavior of the number 



a 2 N (R) = 4<fiS(0)(R/D) 2 + lower order terms, (1) 

where D is an effective length scale, taken here to be the average diameter (D = 7 sm aii-Dsmaii+ 
7iarge-Diarge) of the particles, and R is the radius of a circular observation window. It follows 
from (prj that any point pattern with S(0) = does not possess infinite-wavelength local 
number density fluctuations and thus is hyperuniform. For heterogeneous media consisting 
of finite-volume inclusions, such as the binary MRJ packings we study here, one should 
consider the two-point probability function S^r), defined previously, and the associated 
spectral density x(k) = S^S^O") — <p 2 }(k). Using these functions, the variance o~ 2 (R) in the 
local volume fraction r, the volume fraction of a reference phase in an observation region, 



can be asymptotically written in Mr as [11] : 



a 2 r (R) = [p/(4(j))]x(0)(D/R) 2 + lower order terms. (2) 

As for point patterns, any heterogeneous medium with x(0) = is hyperuniform and lacks 
infinite-wavelength local-volume-fraction fluctuations. Differentiating between the number 
variance and local-volume-fraction fluctuations is essential to characterize MRJ packings 
with a size distribution, and such a distinction has not previously been made in the literature. 

Both the structure factor and spectral density can be numerically obtained by di- 
rect Fourier transform 19] according to S'(k) = A r_1 |^^ 1 exp(ik • r ^ ) | 2 and x(k) = 
V~ 1 \^J =1 exp(i\i-Tj)rh(k', Rj)\ 2 , where rh(k\ R) = ${Q(R— ||r||)}(fc) is the Fourier transform 
of the indicator function for a disk of radius R, Q(x) is the Heaviside step function, and 
{ij} denotes the particle locations. This result can be generalized to include particles of 
any geometry by changing the form of m; the important restriction is that it will only apply 
for media consisting of impenetrable inclusions. The reciprocal vectors k are discretized 
according to the shape of the simulation box. For a square simulation box of side length L, 
the wave vectors are given by k = (27r/L)n for n e 1?. To obtain spherically-symmetric 
forms of the structure factor and spectral density, we radially average over all wavevectors 
within some small spherical shell in reciprocal space. We have considered system sizes up to 
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FIG. 1. (Color online) (a) Structure factor and local configuration (inset) for the binary MR J 
packing. Large and small particles are shown in blue and red, respectively, (b) Corresponding 
spectral density with induced local- volume- fraction fluctuations (inset). The local- volume- fraction 
fluctuations decay logarithmically faster than the volume of an observation window, implying both 
hyperuniformity and the presence of quasi-long-range correlations. 

N = 1000000 particles to probe effectively the small-wavenumber region with even higher 
resolution, due to dimensional scaling, than previous studies jsj. 

Figure [T] depicts both the structure factor and spectral density of the binary MRJ disk 
packing. The small-A; behavior of S(k) is nonvanishing at the origin, showing that the point 
pattern generated by the disk centers is not hyperuniform, a result of the binary nature of 
the packing. Figure [1] also shows that close-packed clusters of small particles permeate the 
system, but these local clusters are separated from each other by effective "grain boundaries" 
arising from the inclusion of large particles. These grain boundaries preclude the suppres- 



d be contrasted with 
5], where S(k) ~ k as 



sion of infinite-wavelength local density fluctuations. This result shou 
previous calculations for MRJ packings of monodisperse spheres in 3D 
k — > 0. It follows that the Torquato-Stillinger conjecture as originally stated for the number 
variance cannot hold for strictly jammed, saturated sphere packings with a size distribution. 
Nevertheless, the physical motivation behind this conjecture is intuitively appealing, and our 
results for the spectral density in Fig. [T]show that it continues to hold in the aforementioned 
generalized sense upon analyzing local-volume-fraction fluctuations. 

We have fit the small-fc region of the spectral density with a minimal third-order poly- 



FIG. 2. (Color online) Jamming constraints on the void space. The left panel shows the MRJ 
state, which is dominated by three-particle loops of the types shown. The right panel is a nearly- 
jammed configuration with volume fraction (f> = 0.75, which exhibits an increasing fraction of 
skewed higher-order loops (also shown) as the void space becomes more connected. 



nomial of the form x(k) ~ X^=o a «^ 1 an< ^ have found a = (1.0 ± 0.2) x 10 5 , strongly 
indicating that this packing is hyperuniform. Thus, like 3D MRJ monodisperse sphere 



packings the spectral density is linear at the origin, implying that asymptotic local- 
volume-fraction fluctuations decay logarithmically faster than the observation window vol- 
ume, i.e., of (-R) ~ (bo ln(i?) + bi)/R 3 as R — > +oo. This behavior implies that S^r) contains 



a QLR power-law tail r 3 , extending to r in R d 11]. 

To clarify the disparity between the number variance and local-volume-fraction fluctu- 
ations, we have examined the distribution of the void space. Since the spectral density 
measures both the inclusion phase and the surrounding matrix jl4j . hyperuniformity is in- 
variant to the choice of the reference phase, and homogeneity in the void phase is sufficient 
to induce hyperuniformity. Figure [2] highlights the void phase for the binary MRJ packing 
in addition to a nonhyperuniform "nearly jammed" disk packing. The MRJ contact net- 
work is dominated by local three-particle contacts (or loops), but by moving away from the 
jammed state, higher-order loops become increasingly common. This observation implies 
that the void-space distribution is inherently constrained by strict jamming. It is true that 
for any strictly jammed, saturated packing of hard disks (2D) the matrix phase must be 
disconnected (see Fig. [21), which, although not essential for hyperuniformity, implies that 



the void distribution is determined completely by the contacts among the particles 24J. In 



a companion paper 



251 ] . we quantify the relationship between jamming and the constrained 




FIG. 3. (Color online) Spectral densities for binary MRJ packings of (left) superdisks (p = 1.5) and 
(right) ellipses (a = b/a = 1.1). Also included are portions of the jammed configurations (insets). 
Large and small particles are in red and blue, respectively. 



void space by measuring the distribution of pore sizes and by utilizing rigorous bounds that 
account for jamming and hyperuniformity. 

Jamming is a critical factor for hyperuniformity to hold; even upon moving infinitesimally 
away from jamming, hyperuniformity is broken. This behavior implies that knowledge of 
local n-particle contacts is equivalent to specifying the distribution of void sizes. Relaxing 
the jamming constraint generates nonuniform, connected void-space regions that eventually 
percolate, skewing the distribution of pore sizes. Thus, the local void sizes in the MRJ 
packing are inherently correlated with each other via the contact network. Specifically, the 
presence of a large void limits the sizes of neighboring voids since the particles are strictly 
jammed. The fact that hyperuniformity is lost upon breaking the contact network implies 
that these correlations are inherently quasi-long-ranged, thereby suggesting a fundamental 
relationship between the distribution of the void space and spatially-extended correlations. 

We have extended this analysis to MRJ packings of noncircular and nonspherical particles. 
Figure |3] shows the spectral densities associated with binary MRJ packings of hard superdisks 



and ellipses [26|. In both cases the packings are hyperuniform even though the inclusion 
shapes are no longer isotropic. Although the local contacts for these packings are less 
uniform than for hard disks, the constraints of saturation and strict jamming again limit 
the distribution of pore sizes and enforce hyperuniformity. Also, the small-fc behaviors of 
the spectral densities for both the superdisk and ellipse packings are linear. 



It has been previously observed that this small-wavenumber scaling is directly related to 
the "degree of order" of the packing |5|, |27|. Linear small-/c behavior in the spectral density 
therefore reflects a reasonable amount of variability in the shapes and sizes of the local 
voids. Lower-order small-wavenumber terms in the spectral density (e.g., x(k) ~ k 1 ^ 2 ) are 
associated with even greater variability in the distribution of the void sizes and therefore 
reflect increasing probabilities of observing large voids. These types of local- volume-fraction 
fluctuations are thus incompatible with strict jamming, and the small-wavenumber scaling of 
the spectral density is minimized to the smallest positive integer value [x{k) ~ k]. Our results 
thereby suggest that QLR correlations in S^r) as induced by a linear small-wavenumber 
region of the spectral density are a likely sufficient indicator of MRJ particle packings. 

We have established universal features of MRJ particle packings using local-volume- 
fraction fluctuations. Our results generalize the QLR pair correlations of 3D MRJ sphere 
packings to saturated, strictly jammed hard-particle packings and provide insight into the 
origin of these correlations. These packings are characterized by QLR pair correlations 
^ r -{d+i) j n j^d^ w j^j 1 a ii near small-wavenumber region in the spectral density, and moving 
even slightly off the jamming point destroys hyperuniformity. We have justified this prop- 
erty using the variability of the void space, which is both constrained by the local contact 
network and maximally disordered among all related saturated, strictly jammed packings. 
Our arguments incorporate higher- dimensional packings, and future work will examine MRJ 



packings of ellipsoids 



281 ] and superballs 29j] . We have also shown that local-number-density 



fluctuations of the particle centers are in general not sufficient to characterize MRJ pack- 
ings because the inclusion shapes and sizes are rigorously linked to the void space, meaning 
that the appropriate quantities to examine are the two-point probability function So(r) 
and the associated spectral density x(k). We mention that recent, independent work [30] 
has examined fluctuation-response relations in MRJ packings of polydisperse spheres with 
similar conclusions to our own. However, no mention is made of quasi-long-range corre- 
lations, and MRJ packings of nonspherical particles are not considered. Our arguments 
provide insight into the structural properties of jammed particle packings and the nature 
of the jamming transition and suggest certain quantum many-body systems and cosmo log- 
ical structures with special QLR correlations iMSJ] are statistically "rigid" with vanishing 
infinite-wavelength number density fluctuations vanish and nonanalyticity in the structure 
factor at small wavenumbers. 
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